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Abstract 


A new numerical scheme based on Genocchi polynomials is constructed 
to solve fractional Sturm—Liouville problems of order two in which the 
fractional derivative is considered in the Caputo sense. First, the differen- 
tial equation with boundary conditions is converted into the corresponding 
integral equation form. Next, the fractional integration and derivation op- 
erational matrices for Genocchi polynomials, are introduced and applied 
for approximating the eigenvalues of the problem. Then, the proposed 
polynomials are applied to approximate the corresponding eigenfunctions. 
Finally, some examples are presented to illustrate the efficiency and accu- 
racy of the numerical method. The results show that the proposed method 
is better than some other approximations involving orthogonal bases. 
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1 Introduction 


The purpose of solving the Sturm—Liouville problem is to find the eigenval- 
ues and eigenfunctions of a set of equations, which have many applications 
in mathematics and physics. The description of the vibrations of a string 
or a quantum mechanical oscillator is modeled as Sturm—Liouville equations 
(2, 3, 9]. For more information about the application of Sturm—Liouville 
problems, we refer to [4, 5, 9, 10, 14] and references therein. With the ad- 
vent of fractional calculus [8, 20, 24], the use of fractional derivatives in 
Sturm-Liouville equations led to a new class of equations, which are known 
as fractional Sturm—Liouville problems. Various types of fractional order 
derivatives are defined [11, 12, 18, 21, 33], but most researchers use Caputo’s 
concept for Sturm—Liouville equations because Caputo fractional derivative 
is more compatible with practical application in physics and engineering 
[22, 23, 25, 26, 30]. The analytical solution to Sturm—Liouville problems 
is usually not computable. This problem has limited the application of these 
equations in various fields. In simpler cases, the analytical solutions to this 
equation can be expressed in terms of specific functions, such as Mittag-Leffler 
functions, which have their own complexities in terms of calculations. So, in 
most cases, scientists seek to find numerical methods to solve Sturm—Liouville 
problems. 
In this study, we consider the following fractional Sturm—Liouville prob- 
lem: 
0 Dey(a) + (Ar(z) — g(z))y(@) =0, O<e<1, (1) 


where g(a) and r(a) are real value continuous functions and r(x) 4 0 for 
x € [0,1], 1 <a < 2, and {D® denotes the fractional Caputo derivative and 
the boundary conditions are as follows: 


ay(0) + by!(0) =0, 
ee tay(l) 0. (2) 


where a, b, c and d are real constants and a? + b? £4 0, c? +d? £0. 

The authors in [2] established sufficient conditions for the existence and 
uniqueness of solutions for various classes of initial and boundary value 
problems for fractional differential equations involving the Caputo fractional 
derivative. Also, the existence and uniqueness of the solution for a fractional 
Sturm-Liouville boundary value problems based on the Banach fixed point 
theorem were proved in [15]. We also refer to [6, 13, 34] for more studies on 
the existence and uniqueness of the solution for boundary value problems of 
types (1) and (2). 

A wide range of numerical methods has been used to solve problems (1) 
and (2) [1, 4, 5, 7, 18, 19]. The Laplace transform method is applied to 
convert (1)—(2) to the equivalent integral equation with a weakly singular 
kernel in [31]. Then, the authors applied a piecewise Lagrange integration 
method to solve the corresponding integral equation numerically. Inspired 
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by their work, we apply the Genocchi polynomials approximation method to 
the problem (1)—(2). 

Genocchi polynomials, which were introduced in [32], are very impor- 
tant and useful polynomials. These polynomials share some great advan- 
tages with Bernoulli and Euler polynomials for approximating an arbitrary 
smooth function [32]. Genocchi polynomials were applied for solving integer- 
order delay differential equations [16], fractional optimal control problems 
[28], and fractional pantograph equation [16]. The numerical solutions ob- 
tained by Genocchi polynomials are comparable or even more accurate com- 
pared to some well-known existing methods. In this study, motivated by 
these advantages, we define and successfully apply the operational matrices 
of Genocchi polynomials to approximate the eigenvalues and the correspond- 
ing eigenfunctions of the Sturm—Liouville problem (1) and (2). First, a new 
approach for calculating the fractional derivative, integration, and product 
operational matrices is introduced, and then the operational matrices are ap- 
plied for problems (1) and (2). The structure of the paper is as follows: In 
Section 2, we recall some definitions and results related to fractional calculus. 
In Section 3, we construct the new numerical method. Some illustrative ex- 
amples are provided in Section 4. Finally, the conclusion is given in Section 5. 


2 Preliminaries 
In this part, the definition of fractional calculus, Genocchi polynomials, and 
their attributes are explained. 


2.1 Fractional calculus 


The Riemann-Liouville fractional integral o/% of order 0 < a < 1 is presented 
with (see [29]) 


Tay Jo (w@— 8)? +u(s)ds, a>, 


( 3 
u(x), a=0. Y) 


olfu(r) = 


One of the fundamental attributes of the operator 91% is 


The Caputo fractional derivative of order a > 0 is determined as [29] 
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a” 1 = 
Cpa —, pr-a _ n-—a-—1,(n) 
§ Doula) =o Be-° Tule) = ay f (@—sye tw (sds, 6) 


where n is an integer (n—1<a<n) and u™ € L'(0, 1]. 


The main relationship between the Riemann—Liouville integral operator 
and Caputo fractional derivative is as follows: 


0 De off u(x) = u(z), 


plo 2 De al )-5 ul) (0%) (n—1l<a<n), (8) 


where u‘") € E*/0,1], 7 =0,1,...,n—1. 


2.2 Genocchi polynomials and properties 


The Genocchi numbers g,, and the Genocchi polynomials G,,(x) are defined 
as the coefficients of the exponential generating functions as follows: 


2t ee 
—— = nz) t : tr 
ss i “ th<n (7 
= = Gul e. jae (3) 
TT. 
ef+1 = nl’ 
The Genocchi polynomial G,,(x) is a polynomial given by 


n , . . 
5 ah ee ‘y= (jon WS My Dy exo (9) 
k=0 


The Genocchi polynomials have interesting properties, some of which are as 


follows: 
: 2(—1)"n!m! 


: Gy (2)Gal tax = ee a n,m > 1, (10) 
£ Galo) = 0G, 1{2), n= 1, (11) 
ieee, hei. (12) 


The set of Y = {Gi(x), Go(x),...,G@n(x)} C L?[0,1] is a linearly indepen- 
dent set (see[16]). Any f(x) € L?[0,1] has a unique best approximation in 
Span(Y), as f(x), which can be represented by Genocchi polynomials as 
follows: 
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N 
f@) 2 fylz)= >) Gale) =O’ Gia), (13) 

n=1 
where C = [c1,¢2, ...,cw]” is the coefficient vector and G = [Gi(x), Go(z), 


..;Gn(x)|" is the Genocchi polynomials vector. The property of the best 
approximation requires 


and then 
CT = (f(x),G(x))? Wt, (14) 


where 


(f(x), G(x)) = [(f(w), Gi(x))  (f(@), a(x)... (F(2), Gw(a))), 


and W = (G(x), G(x)) isan N x N symmetric matrix, and by (10) its entries 
are calculated as follows: 


1 
2(—1)"n!m! 
= —————gmyn-. (15 


The following lemma provides the upper bound for the error of function 
approximation by Genocchi polynomials. 


Lemma 1. Suppose that f(x) ¢ CN*+![0,1] is approximated by truncated 
Genocchi polynomials f(a) in (13). Then 


R 


If(@) — vl $ Gay ERS (16) 


where R = maxzejo,1 |f%*) (2)|. 


Proof. We refer the reader to [17]. 


2.3 Genocchi operational matrices 


By (11), the derivative of G(#) can be expressed as follows: 


<G(a) = DG(z), (17) 


where 
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000... 0 O00 
200... 0 O00 
030... 0 O00 
p—|004... 0 00 (18) 


000...N—100 
000... 0 NO 


is called the operational matrix of integration for Genocchi polynomials. It 
is obvious that 


—G(«) = D"G(z). (19) 


The operational matrix for fractional order derivative of order a in Caputo 
sense is defined as follows: 


¢ D2G(a) =D" Gia), (20) 
where D® is an N x N matrix. By using (14), we can write 
D* =Q°w', (21) 


where Q® = (D°G(x),G(a)) is an N x N matrix whose entries are defined 
by 


iy = (D°G: (x), G;(x)) 


r=0 
. 3 Verh! | oktr—ade (22) 
— T(k+1-a) Jo 
se i J f| 
>>> TR Yr i>n. 
eS Th leh prt le) — 


It is obvious that the first n — 1 rows for Q® are equal to zero, where n—1 < 
a<n. 


The operational matrix for fractional order Riemann—Liouville integration 
of order a, (0 < a < 1) is defined as follows: 


I Gta) = P"Gla), (23) 
where P® is an N x N matrix. By using (14), we can write 
P* =v°-w"!, (24) 


where V® = (I°G(x),G(x)) is an N x N matrix whose entries are defined by 
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Vij = ("Gi (2), G5 (2)) 


(25) 


_ 3 Tie! 
T(k+1+aj(k+r+1+a) 


The product of two Genocchi vectors can be expressed by Genocchi vector 
as follows: 
G(x)G? (x)C = CG(z), (26) 


where C is called the operational matrix for the product. By using (14), the 
product operational matrix C can be defined as follows: 


c=sw-t, (27) 


where § = (G(x)G7(x)C,G(z)). The entries of § are determined as 


~ x (28) 
=( Giz) >) Gn (2):G;(@) ) = >> cede, 
k=0 k=1 
where 
1 
dijk = Gi(x)G(x)G; (a)da 
0 
: es : (29) 
= by S- a Vr Vs Vt 
r=0 s=0 t=0 r+ts+t+l 


3 Method of solution 


In this section, we apply Genocchi polynomials to solve problem (1). Consider 
the fractional Sturm—Liouville problem (1) with the boundary conditions 
(2). The following Lemma converts this problem to the equivalent integral 
equation. 


Lemma 2. The function y() is a solution for the fractional Sturm-—Liouville 
problem (1)—(2) if and only if y(«) satisfies the following fractional integral 
equation: 
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y(x) = g(a) [ha oly f(@)lea1 + he D oly f(@)lea)— ofr f(@), (30) 
where f(x) = (Ar(x) — q(a))y(a) and 


off f(t) = 7 “(e— tf (t)at. 


The operator D in (30) is the classical derivative operator and g(x), h,, and 
ha are defined considering the values of a, b, c, and d as follows: 


acx — bc d 
= ——— , h, = 1, ho = if 
g(@) ad — be +ac’ owe, ie ere 
adx — bd c 
= ————,, bk =5, ho =1 if 
g(x) = beak qo he =), ifa 40, d40, 
—adx +b 
4 ne, et ifb40, d£0, 
—ad + be — ac d 
—acx + be d : 
ae aa hy =1, ha=- ifb #0, c#0. (31) 


Proof. We refer to [31]. 


To solve the Sturm-—Liouville problem (1)-(2), we apply the Genocchi 
polynomials method to the integral equation (30). First, we approximate 
y(x), r(x), and q(x) in (30) by the truncated Genocchi polynomials (13) as 
follows: 


y(z)~C7G(z), r(z)=RTG(e), g(z)~ Q*G@). (32) 
To apply (32) in (30), we first approximate the elements in (30) as follows: 


f(x) = y(a)(Ar(x) — a(x) 
~ CP G(2)(AG" (z)R — G7 (2)Q) (33) 
= ct QR es Q) G(x) = CTFC(a), 
where R and Q are the product operational matrices corresponding to the 


vectors R and Q, respectively, and F = AR — Q. By (33), the fractional 
integration part of (30) can be written as follows: 


ol f(x) ~ CTFoI2G(x) = CTFP*G(cx), (34) 


where P® is the operational matrix for fractional integration of order a. The 
other part of (30) can be represented as 


D ol f(x) ~ CTFP°DG(x) = CTFP°DG(z), (35) 
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where D is the operational matrix for regular derivative of order one. 


Substituting (32), (33), (34), and (35) in (30) and neglecting the trunca- 
tion errors yield 


CTG(«) = CT { [mF + hoFP°D] G(1)GT — Fp°} G(x). (36) 
The property of linear independency for Genocchi polynomials requires 
(I— A(A))C = 0, (37) 


where 7 - 7 
[A(A)|2 = [miFP° + hoFP°D| G(1)GT — FP*, (38) 


and I is the N x N identity matrix. 


To have nontrivial solutions for (1) and (2) we have to find nonzero solu- 
tions for (37). Therefor, we need to solve the following root-finder problem: 


det (I— A(X)) =0, (39) 


which can be done by mathematical softwares such as MATLAB and Maple. 


After computing eigenvalues, we approximate the eigenfunction y;(t¢) cor- 
responding to the eigenvalue A; by solving the following linear system: 


(I— A(A;))C = 0. (40) 
Since det(I — A(\;)) = 0 and , is the simple root, then we set Cy = c 


and solve (40) to find Cy, k = 1,2,...,N —1 with respect to c. Then the 
eigenfunction y;(t) is obtained by 


HEY CoaGe@,. 40 dyare: (41) 


By using an appropriate auxiliary condition (for example, y;(0) = 1 or 
yi’ (0) = 1), y(t) can be determined uniquely. 


4 Numerical results 


In this section, we present some illustrative examples and show the efficiency 
of the proposed method. 


Example 1. Consider the following fractional eigenvalue problem: 


0 D3 (x) + Ay(z) = O0<e<l, 


(0) = re “0, 2) 
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exact value N = 25 error N = 800 error 
Ai 9.86960440108936 9.869604401089359 2.58e-28 9.86960440 1.09e-9 
A2 39.4784176043574 39.47841760435743 5.22e-19 39.47841760 4.36e-9 
A3 88.8264396098042 88.82643960980423 1.62e-15 88.82643963 2.02e-8 
v4 157.913670417430 157.9136704174555 2.58e-11 157.91367057 1.53e-7 
A5 246.740110027234 246.7401100283036 1.07e-9 246.74011064 6.13¢-7 
re 355.305758439217 355.3057578574429 5.82e-7 355.30575850 6.08e-8 
Az 483.610615653379 483.6106104947717 5.16e-6 483.61061573 7.66e-8 
As 631.654681669719 631.6551696075792 4.88e-4 631.65468191 2.40¢e-7 

Table 2: Eigenvalues for different values of a, Example 1 ( N = 20). 


a= 1.95 
Ai: 9.66077186263 
A2 36.4045254812 
A3  — 80.3290382609 
Ag —-140.076253643 
As: 216.597587060 
As _ 308.346500930 


a=1.9 
9.51414296571 
33.5956740224 
73.0390014094 
124.418589764 
191.145884721 
267.943639675 


a= 1.85 
9.44013733938 
30.9825046888 
66.9448021792 
110.360397400 
170.311838185 
232.015405069 


a=1.8 
9.45685703307 
28.4768769642 
62.2003971381 
97.0631636666 
155.450547906 
196.595011516 


This corresponds to the case a # 0 and c ¥ 0 in (31). Then, we set 
a=c=1,b=d=0, hi = 1, ho = 0, and g(x) = x. For a = 2, (42) 
is converted to the classical Sturm—Liouville problem corresponding to the 
eigenvalues A, = (n7)?. 

In this example, we use the auxiliary condition y/(0) = 1 to approximate 
the eigenfunctions. In Table 1, the numerical results are compared with the 
results of [31]. Table 2 represents the approximated eigenvalues obtained 
by this method. The results in Tables 1 and 2 clearly show that the pro- 
posed algorithm is accurate. Figure 1 shows the first four eigenfunctions for 
a = 1.85, 1.9, 1.95, 2 and N = 15. In Figure 2, the graphs of same order 
eigenvalues are compared for different values of a. 


Example 2. For the second example, we solve the following eigenvalue prob- 
lem: 


A 


Cpa 
D = 
§ Deva) + Goo 


O0<a<l, 
(43) 


For a = 2, the eigenvalues of (43) are An = (nm/In2)? + 1/4, and the 


nm In(1+2) 
In 2 


Cn is constant coefficient [31]. The entries are chosen hy = 1, hg = 0, and 
g(x) = a. We use the auxiliary condition y/(0) = 1 to compute four first 
eigenfunctions of Example 2, which are plotted in Figure 3. Tables 3 and 4 
represent the numerical results obtained for Example 2, which were compared 
with the results in [31]. 


corresponding eigenfunctions are yp,(@) = CnV1+ sin ( ); where 
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a=2 a= 1.95 
0.34 Lgl 
0.24 024 
0.14 o14 
o4 
o-4 
-0.14 
-o.14 
0 0.2 04 0.6 08 1 0 0.2 04 0.6 08 1 
x # 
—h’—1, 4, ——Ay —h,—4, 4, —h, 
a=19 = 1.85 
0.25 
0.24 0.20 
O15 
O14 0.10 
0.05 
o4 
0 
-0.05 
-0.1+ 
0 02 04 0.6 08 1 0 0.2 04 0.6 08 1 
x x 
—’A—h, 4, ——Ay —A, —aA ds; —h, 


a=18 
0.25 
0.20: 
0.15 
0.10: 
0.05 
0 
0.05 

0 0.2 04 0.6 0.8 1 
x 
—}, 24, — 4, —), 


Figure 1: First four eigenfunctions for different values of a, Example 1. 


Example 3. For the next example, we solve the following eigenvalue prob- 
lem: 
6 Dey(x) + (A+ 10sin rx)y(ax) = 0, 0<2r<l, 


(O)=0, 4-0 4) 
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0.15 


132 


04 


0.6 


0.6 08 
x E 
o=2 o=1.95 o=1.9 o=1.85 oO=1.8 o=2 a= 1.95 o=1.9 O=1.85 O=1.8 
as 
0.10: 
0.05: 
0 
0.05, 
0.10: T T T T 
0 0.2 04 0.6 0.8 1 
‘ 
o=2 v1.95 O=1,9 O=1,85 O18 o=2 o=1.95 o=1.9 O=1,85 O18 
Figure 2: First four eigenfunctions for Example 1. 


Table 3: Eigenvalues for a = 2, Example 2. 


proposed method 


method of [31] 


exact value 


N = 25 


error 


N = 500 


error 


AL 20.7922884552238 
dA2 82.4191538208953 
A3 185.130596097014 
da 328.926615283581 
As 513.807211380596 
ro 739.772384388058 
Az 1006.82213430597 
As8 1314.95646113432 
Ag 1664.17536487313 
A1o 2054.47884552238 


20.7922884552238 
82.4191538208952 
185.130596097016 
328.926615287143 
513.807211123754 
739.772393148381 
1006.82198782297 
1314.95690886405 
1664.19294987497 
2054.28429017518 


2.70e-17 

6.85e-14 

1.91e-12 
3.56¢e-9 
2.57e-7 
8.76¢e-6 
1.46e-4 
4.48¢e-4 
1.76e-2 
1.95e-1 


20.79228809 
82.41914815 
185.13056780 
328.92652741 
513.80700157 
739.77196118 
1006.82137649 
1314.95521911 
1664.17346838 
2054.47611287 


3.65e-7 
5.67e-6 
2.83e-5 
8.79e-5 
2.10e-4 
4.23e-4 
7.58e-4 
1.24e-3 
1.90e-3 
2.73e-3 


Setting hy = 1, hg = 0, and g(x) = x, we use the proposed method to 
solve (44). With N = 20 and 40 digits for computation in Maple software, 
some first eigenvalues are listed in Table 5 for different values of a. 

We use the auxiliary condition y;(0) = 1 to compute four first eigen- 
functions of Example 3, which are plotted in Figure 4 for different values of 


Qa. 


Example 4. Consider the following eigenvalue problem: 
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Table 4: Eigenvalues for different values of a, N = 20, Example 2. 


a= 1.95 


a=1.9 


a= 1.85 


a=1.8 


20.6677971198 
75.7710019802 
167.652678491 
291.242134899 
451.065004856 
641.015025258 
867.315229409 
1122.28888722 
1414.21783814 
1732.91424589 


20.7298859439 
69.4887851211 
153.001698048 
257.628271127 
398.913421123 
554.923522781 
752.921297326 
956.431286954 
1212.08349835 
1458.60978678 


21.0420113791 
63.3252982915 
141.496620930 
226.249088945 
358.715501294 
474.802195243 
669.708580955 
798.315596410 
1084.45731615 
1178.69721277 


21.7289820388 
56.9307959084 
134.502772594 
193.872634358 
345.306193742 
379.243396666 
2260.68651346 
3571.46938188 
8833.11202980 
80410.2798175 


x 


a=2 


o=1.95 elo 


O=1.85 o=1.8 


e 


02 04 


o=2 OF1.95 


= 1.85 O18 


a=2 


o=1.9 


o=1,.85 


0.2 04 


a=2 o=1.95 


O=1,85 o=1.8 


Figure 3: First four eigenfunctions for Example 2. 


6 Doy(x) + (2re” — 5sin ra) y(x) = 0, 
y(0) = y'(0) £0, y(1) =0. 


O0<a<l, 


(45) 


We set a = 1, b= -1,c = 1, andd=0. Then hy = 1, ho = 0 and 
(x + 1)/2 (see (33) ). With N = 25 and 40 digits for computation, 


some first eigenvalues are presented and compared in Table 6 for different 
values of a. 
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Table 5: Eigenvalues of Example 3 (N = 20). 
a=2 a= 1.95 a=1.9 a=1.8 
AL 18.3200943208 18.1463655269 18.0384472376 18.0792788380 
A2 46.2367305574 43.1880344976 40.4110340394 35.3810395900 
AZ 95.3857838965 86.8856296038 79.5840346405 68.6611728543 
Aa 164.391065891 146.560791927 130.915815214 103.643270952 
A5 253.178225348 223.035468745 197.577387685 161.771332663 
6 361.722182212 314.767445388 274.375010614 203.160705578 
Az 490.014088674 423.155049789 367.477687375 305.476138569 
Ag 637.940873962 546.469154163 468.674770446 315.791118150 
A a, 
0.15 
0.10: 
0.05, 
0 
0.05 
0.10 
0 -0.15 
Oo 0.2 04 06 O08 1 0 0.2 04 0.6 08 A 
a= = 1.95 o=1.9 = 1.85 o=1.8 a2 0=1,95 o=1.9 O=1.85 a=18 
a, A, 
0.10: 
0.06 
0.05 0.04. 
0.0: 
0 0 
-0.0; 
0.05 0.04. 
0.06 
0.10, 
0 0.2 04 0.6 08 1 0 02 04 0.6 08 1 
a2 o=1,95 o=1.9 o=1.85 o=1.8 o=2 o=1.95 o=1,9 o=1,85 o=1.8 
Figure 4: First four eigenfunctions for Example 3. 
sis ae P _ P = ; 
We use the auxiliary condition yj(1) = 1 and yj(1) = —-1, 7 # 1 to 


compute four first eigenfunctions of Example 4, which are plotted in Figure 


5. 


Example 5. Consider the following forth-order eigenvalue problem: [27] 


v (t) — (0.02t?0' (t))’ + (0.0001t* — 0.02)u(t) = A v(t), 


v(0) v(1) = v"(1) 


=v"(0) =0, 


= 0. 


ves, 


(46) 


The eigenvalues of this problem are the square of the eigenvalues of the 
following second order Sturm-Liouville problem (see [27]): 
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Table 6: Eigenvalues of Example 4. 


proposed method, N = 25 


a= 2 a= 1.95 a=1.9 a= 1.85 a=1.8 
Ai 2.6213001864 2.5785602459 2.5408885358 2.5083120348 2.4809701676 
dA2 8.3495214073 7.7826681845 7.2770461324 6.8263714594 6.4250121537 
A3 20.160277342 18.278658822 16.631507344 15.191185649 13.934739577 
Aa 37.779499263 33.634640153 30.049170175 26.945962302 24.258034906 
A5 61.245497320 53.803913449 47.432736730 41.977052803 37.309595288 
6 90.568134413 78.731592840 68.678948337 60.132525229 52.854782758 
Az 125.75080891 108.37362381 93.720958762 81.357289626 70.931030267 
Ag 166.79492842 142.69052813 122.48797531 105.53364534 91.272428145 
method of [31], N = 800 
a=2 a= 1.95 a= 1.85 
Ai 2.62130018 2.57856033 2.50831250 
dA2 8.34952147 7.78266742 6.82636707 
A3 20.16027739 18.27866231 15.19120752 
ra 37.77949944 33.63463032 26.94588552 
As 61.24549782 53.80393654 41.97727381 
Ay a, 
0.44 
0.44 0.34 
0.24 
0.34 
0.14 
0.27 
o 
0.14 
-0.174 
0 -0.24 
0 0.2 0.4 0.6 0.8 1 0 0.2 04 0.6 0.8 1 
o=2 O=1.95 o=1.9 o=1.85 o=1.8 a=2 o=1.95 o=1.9 o=1.85 o=1.8 
ay dg 
0.24 0.54 
0.44 
0.14 
0.34 
i 
0.24 
-0.14 0.14 
-0.24 07 
on 
0.34 
0.24 
0.44 -034 
oO 02 04 06 08 i 0 02 04 016 Og H 
o2 O=1.95 o=1.9 o=1.85 o=1.8 o2 o=1.95 o=1.9 o=1.85 o=1.8 
Figure 5: First four eigenfunctions for Example 4. 
—u" (t) + 0.01t?u(t) = » u(t), 0=¢<5 (47) 


By taking the variable change t = 5a, the following problem is achieved: 
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y” (x) +(25d — 6.25a7)y(x) = 0, 
y(0) =0, y(1) =0. 


Taking hi = 1, he = 0, g(x) = 2, r(x) = 25, and q(x) = 6.2527, with N = 
25 and 40 digits for computation in Maple software, some first eigenvalues 
are presented and compared with the results of [27] in Table 7. We use the 
auxiliary condition y/(0) = 1 to compute four first eigenfunctions of Example 


(eee ii, 
aaa (48) 


5, which are plotted in Figure 6. 


Table 7: Eigenvalues of Example 5. 


proposed method 


[27] 


AL 0.215050864368 0.2150508644 
A2 2.754809934683 2.7548099347 
A3 13.21535154056 13.2153515406 
4 40.95081975918 40.9508197592 
A5 99.05347806596 99.0534780635 
ro 204.3557315637 204.355732268 
Az 377.4304091791 377.430420689 
As 642.5918663327 642.590868170 


5 Conclusion 


In this paper, the eigenvalues of the second-order fractional Sturm—Liouville 
problem were approximated by using Genocchi polynomials. The opera- 
tional matrix for fractional integration, fractional derivative, and product 
was evaluated. Then, the operational matrices were applied to the fractional 
Sturm-Liouville problem to convert it into a homogeneous system of lin- 
ear equations. The eigenvalues of the coefficient matrix corresponded to the 
eigenvalues of the main Sturm—Liouville problem. The presented illustrative 
examples verified that the proposed method generated more accurate approx- 
imations compared to the results, which were reported in other papers. Like 
other methods in the literature, the results were accurate for lower indices 
but not very accurate for higher indices. After computing the eigenvalues, 
the corresponding eigenfunctions were approximated by applying an auxil- 
iary condition. The results were in good agreement with the exact solutions 
or the results of other papers. 
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